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1.  Introduction 


The  von  Kirm&n  equations  for  a  damped  plate  are  given  by  (see,  e.g.,  [7]  or  [9]) 

Aalh  +  ^[lk>V>a]  =  0  in  ft 

and 

A2V>a  -  =  A$  in  ft 

where 

,  m  -  _  a2v>  d*<t> 

6x\dx\  8x\dx\  dllX2  &ZlX2 

Here,  ft  is  a  bounded,  convex  polygonal  domain  in  R,2,  in  denotes  the  Airy  stress  function,  V’z 
denotes  the  deflection  of  the  plate  in  the  direction  normal  to  the  plate,  and  A g  is  an  external  load 
normal  to  the  plate  which  depends  on  the  loading  parameter  A. 

The  boundary  conditions  on  T  =  dCl  are  given  as  follow: 

*  =  f!r='k  =  i£=0  onr’ 

where  9(-)/dn  denotes  the  normal  derivative  in  the  direction  of  the  outer  normal  to  I\ 

For  reasons  to  be  explained  later  we  introduce  appropriate  rescalings,  i.e.,  by  replacing  in  by 
Ay»i,  i> 2  by  A xfn,  and  g  by  A g.  Then  we  can  rewrite  the  Kirm&n  equations  as  follows: 


A2lh  +  =  0  in  ft. 

(1.1) 

A3tfe  -  A[t&i,^a]  =  A$  in  ft, 

(1.2) 

Vl  =  ^='6j  =  ^=0  °"r- 

(1.3) 

We  introduce  the  standard  Sobolev  spaces  (see,  e.g.,  [1]) 

#2(ft)  =  e  ff2(ft)  |  i>  =  o,  =  o  on  r}  , 

H02(ft)=(H0a(ft))2,  H->(Sl)=(H$(Sl))\  and  ff-2(ft)  =  (H02(ft))*, 

the  bilinear  form 

<*(tM)=  /  AV>A<Mft  Vi>,<f>€  H2(H) 

Jo 

and  the  L2(ft)-inner  product 


(V>,^)=  /  i><t>d&  i  i>,<f>€  Zr3(ft) . 

Jo 

Then  we  may  define  the  following  weak  formulation  of  the  von  Karmin  equations  (1.1)-(1.3):  And 
+  =  (th ,  ih)T  €  H*(Sl)  such  that 

a(in,<h)+^([i>2,H,<tn)=0  V<tn€H*(Sl)  (1.4) 
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and 

«(V*,fc)-A([V»i,M*)  =  Ms,<h)  V<h£H$(Sl).  (1.5) 

Here  the  notation  (•,  •)  stands  for  the  L2(fl)-inner  product  and  (•,  •)  the  duality  pairing.  Using  the 
identity 

(hM.C)  =  (ty,C U)  VtM,C€  (1.6) 

one  can  show  that  for  each  g  6  /?-2(Q),  (1.4)-(1.5)  possesses  at  least  one  solution  ^  =  (0i,  *h)T  € 
Hq(U)  and  that  all  solutions  of  (1.4)-(1.5)  satisfy  the  a  priori  estimate 

11^113  + l!V>3||2<C||S||_2;  (1.7) 

see,  e.g.,  [9],  for  details.  In  the  sequel  a  solution  to  (1.1)-(1.3)  will  be  understood  in  the  sense  of 

0-4M1-5). 

Given  a  desired  state  =  (V'io,  V'ao)  6  L2(fi),  we  define  the  functional 

JO>,9)  = 

fn ( (( 'h  ~  ^1°)J  +  (to  ~  too)2)  dSl  +  ^J^g2<Kl 

for  all  r/>  =  {1n,>h)T  €  and  g  €  L2(f l).  We  then  consider  the  following  optimal  control 

problem  associated  with  the  von  Karman  plate  equations: 

min{J(1>,g)  |  g  €  0}  subject  to  (1.4)-(1.5),  (1.9) 

where  0  is  a  convex,  closed  subset  of  L2(Sl).  The  physical  interpretation  of  this  optimal  control 
problem  is  that  we  wish  to  match  a  desired  Airy  stress  distribution  and  a  deflection  distribution 
on  the  entire  plate  by  choosing  an  (optimal)  external  load  g  from  an  admissible  set  6. 

Our  plan  is  as  follows.  In  §2,  we  show  that  an  optimal  solution  exists.  In  §3,  we  use  Lagrange 
multiplier  techniques  to  derive  an  optimality  system  from  which  optimal  states  and  controls  can  be 
deduced.  We  then  specialize  to  the  case  0  =  L2(0)  and  derive  an  optimality  system  of  equations. 
In  §4,  we  define  conforming  finite  element  approximations  of  the  optimality  system  of  equations 
and  derive  error  estimates. 

2.  Existence  of  an  optimal  solution 

We  will  prove  in  this  section  that  an  optimal  solution  exsits  for  the  minimization  problem 

(19). 

Theorem  2.1.  There  exists  a  (<f>,g)  6  x  0  such  that  (1.8)  is  minimized  subject  to 

(1.4M1.5). 

Proof.  We  choose  a  g(°)  €  0  and  let  V’2°^)T  €  Hq(SI)  be  a  solution  of 

aU>[°\<h)  +  =0  V  4>\  €  ff2(0) 

and 

«(V40),  to)  -  A(tyi0),  40)],  to)  =  a (s<0),  <h)  v  <h  €  . 
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The  existence  of  such  a  was  established  in  [9].  We  see  that  (V^0* ,  i^0))  €  x  0  satisfies  the 

constraints  (1.4)-(1.5)  and  J <  oo.  We  also  note  that  J (f,g)  >  0  for  all  (f,g)  e  XxQ. 
Thus  we  may  choose  a  minimizing  sequence  {(^"* ,  p(n))}  C  Ho(ty  x  ®  such  that  for  some  constant 
M  >0, 

J{+(n\g^)<M,  (2.1) 


and 


a(Wn\<M  +  =  0  V<fc  €  H2(Sl) 


a(1>in),<h)  -  A([*Sw\*$n)],fc)  =  \{g<n\<h)  V<fc  G 


(2.2) 


(2.3) 


Using  (2.1)  and  the  definition  of  the  functional  J  we  deduce  that  the  sequence  {|)y*n)||o}  is  bounded 
so  that  by  the  apriori  estimate  (1.7),  the  sequence  {HV^Ih}  >s  bounded  as  well.  Hence  by  choosing 
a  subsequence,  we  have  that 

g(n)  g  in  £,2(11) 


and 


in  H$(Sl) 


for  some  g  €  L2(Sl)  and  iff  €  Hq(SI).  Using  the  convexity  and  closedness  of  9,  we  deduce  that 
j  €  0.  By  the  compact  imbedding  Hq(SI)  L2(H),  we  obtain  — ►  iff  in  L2(Sl)  so  that  for 

each  ^  =  (<f>i,<h)T  G  Hq(H),  we  have  that  (also  using  the  continuous  imbedding  Hq(U)  «-*  L°°(Sl)) 


=  .'22, 5 


=  ^(ltk><fo]>tfe)  =  ^([V>2,tfo],^l) 


and 


Jim  -  ([i>[n),^{2n)},<h)  =  Jim  -  ([^{n\<fe],v4n)) 

Hence  we  may  pass  to  the  limit  in  (2.2)-(2.3)  to  show  that  (^,p)  satisfies  (1.4)-(1.5). 

Finally,  we  use  the  sequential  lower  semi-continuity  of  the  functional  J  to  obtain  that 

J(+,g)<  lim 

n— *oo 

Hence,  we  have  proved  that  a  solution  (1>,g)  G  U<j(H)  x  9  exists  that  minimizes  (1.8)  subject  to 
(1.4)-(1.5).  0 


3.  An  optimality  system  of  equations 

In  this  section  we  assume  that  (1>,g)  G  #o(fi)  x  9  is  an  optimal  solution  for  the  minimization 
problem  (1.9)  and  we  attemp  to  characterize  the  optimal  solution  as  the  solution  for  a  system 
of  partial  differential  equations.  To  be  precise,  we  use  Lagrange  multiplier  rules  to  derive  an 
optimality  system  of  equations. 
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We  define  the  Lagrangian  for  the  constrained  minimization  problem  (1.9)  as  follows: 


£(*.«.■?)=  { a(V»i , *71 )  + 

+  o(V>2, rja)  —  X([iln,ih],V»)  - 

for  all  (*,S,*)  €  N0a(ft)  x  L>(Sl)  x  ff02(«)- 

By  formally  taking  variations  in  the  Lagrangian  with  respect  to  ip  and  g ,  we  obtain: 

«(Ci,»h)  -  A([ik,t»],Ci)  =  Kin  -  V>io,Ci)  VCi  € 

«(Ca , »?2 )  4-  A([^,tji],<2)  -  A([V»i , *72], C2) 

=  A(V>j-V>20,C2)  VC2  €H02(ft) 


and 


^  (z,  z)  +  A(z,  Tfi)  -  |  (5,  g)  -  X(g,  qj)  >  0  V  z  €  0 . 


(3.1) 

(3.2) 

(3.3) 


For  each  e  €  (0, 1)  and  each  t  €  9,  we  set  z  =  ct  +  (1  -  t)g  £  8  in  (3.3)  to  obtain 
^  (t  -  g,t  -  g)  +  c(t  -  gtg)  +  e(t  -  g,rn)  >0  Vt  €  6 

so  that,  after  dividing  by  e  >  0  and  then  letting  (  -*•  0+,  we  obtain 

(t-9,9  +  *»)  >  0  vtee. 


(3.4) 


Now  we  show  that  there  does  exist  an  9  (SI)  satisfying  (3.1)-(3.3),  or  equivalently,  (3.1)- 
(3.2)  and  (3.4),  so  that  we  are  justified  to  compute  a  triplet  (ip,g,r))  from  (1.4)-(1.5),  (3.1)-(3.2) 
and  (3.4).  In  this  paper  we  will  not  address  the  uniqueness  of  solutions  for  the  system  (1.4)-(1.5), 
(3.1>(3.2)  and  (3.4). 

We  first  quote  the  following  abstract  Lagrange  multiplier  rule  whose  proof  can  be  found  in 

IW- 

Theorem  3.1.  Let  Xi  and  Xj  be  two  Banach  spaces  and  &  an  arbitrary  set.  Suppose  J  is  a 
functional  on  X 1  x  0  and  M  a  mapping  from  Xj  x  0  to  Xj.  Assume  that  (u,g)  is  a  solution  to 
the  following  constrained  minimization  problem: 


M(u,g)  =  0  and  there  exists  an  «  >  0  such  that  J(u,g)  <  J(v,z ) 
for  all  («,z)  such  that  J|u  -  t?|jxj  <  <  and  M(v,z )  =  0. 


(3.5) 


Let  U  be  an  open  neighborhood  of  u  in  X\.  Assume  further  that  the  following  conditions  are 
satisfied: 

(A)  for  each  z  6  0,  t)t-*J(t),  z)  and  v  M(v,z)  are  Prechet-differentiable  at  v  =  u ; 

(B)  for  any  v  €  U,  zj,  Z2  €  9,  and  7  €  [0, 1],  there  exists  a  z7  =  z~,(y,  zi,  zi)  such  that 

M(v,Zy)  =  7Af(t\zi)  +  (1  -  7 )Af(v,zi) 
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and 

(C)  Range(M*(u,  g))  is  closed  with  a  finite  codimension\ 
and 

(D)  the  algeraic  sum  Mu(u,g)Xi  +  M(tt,  0)  contains  0  6  Xj  as  an  interior  point, 

where  Mu(a,g)  denotes  the  Frechet  derivative  of  M  vrith  respect  to  u.  Then,  there  exists  a  a  ij  €  Xf 
such  that 

( r h  Mu(u,  g)v)  -  < Ju(«,  g),  t>)  =  0  V  v  €  X\ 

(or  equivelently,  [Mtt(ii,  p)]*q  -  J„(it,$)  =  0)  and 

min  jC( u,  z ,  tj)  =  £(«,  g,  q)  , 

where  C(u,g,Tf)  =  J(u,g)-(rj,  M(u,g))  is  the  Lagrangian  for  the  constrained  minimization  problem 

(2.5)  and  where  Ju(a,g)  denotes  the  Frechet  derivative  of  J  with  respect  to  u.  □ 

We  now  define  some  spaces  and  operators  in  order  to  rewrite  the  constraint  equations  (1.4)- 

(1.5)  into  a  form  that  will  facilitate  our  verification  of  the  assumptions  in  Theorem  3.1. 

We  define  the  spaces  X  =  Y  =  G  =  L2(n),  and  Z  =  L1(ft).  By  compact 

imbedding  results,  Z  <-*<-+  Y .  The  control  set  0  is  a  closed,  convex  subset  of  G  =  L2(ft).  We  also 
note  that  Y  =  X*. 

We  define  the  continuous  linear  operator  A  €  £(X;  Y)  as  follows:  for  ^  €  X  =  H${ ft), 
A+  =  /  €  Y  =  ff-2(ft)  if  and  only  if 

*(il>u<h)  =  (fi,<h)  V^i  eff0J(fi) 


and 

a(ih,<h)  =  (fi,<h)  6  Bq(U)  . 

It  can  be  easily  verified  that  A  is  self-adjoint,  invertible,  and  A~l  €  C(Y ;  X). 
We  define  the  (differentiable)  nonlinear  mapping  N  :  X  — ►  V  by 


or  equivalently 

<*(*),  *)  =  ^  (fa ,  M  )  -  (fa ,  M  4a)  v +  =  (<h,<h)T  ex 

and  define  K  :  g  €  L2(fl)  -♦  Y  by 

*-(!)• 

or  equivalently, 

(tf5,4)  = -<sfa>  *4  =  fafa)€X. 
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Clearly,  the  constraint  equations  (1.4)-(1.5)  can  be  expressed  as 

A++\N{*)  +  \Kg  =  0. 

We  are  now  in  a  position  to  prove  the  existence  of  a  Lagrange  multiplier  q  €  that 

satisfies  (3.1)-(3.2)  and  (3.4). 

THEOREM  3.2.  Assume  (^,$)  €  H$(Sl)xQ  is  an  optimal  solution  that  minimizes  (1.8)  subject 
to  (1.4)-(  1.5).  Then,  for  almost  all  A  €  A,  there  exists  a  Lagrange  multiplier  q  €  /7<j(0)  satisfying 
the  Euler  equations  (3.1)-(3.2)  and  (3.4). 

Proof.  We  define  a  mapping  M  :  Hq(SI)  x  9  — ►  H~2(Sl)  (i.e.,  X  x  9  — +  V)  by 

M(*,  z)  =  A+  +  A  N{4>)  +  A  Kz  V  (*,  z)  e  H*({1)  x  0 . 

We  see  that  minimizing  (1.8)  subject  to  (1.4)-(1.5)  can  be  stated  as  minimizing  J{1>,g)  subject  to 
M(4,g)  =  0.  We  first  verify  that  the  hypotheses  (A)-(C)  of  Theorem  3.1  hold  with  Xi  =  X  and 
Xi  =  Y. 

Obviously,  for  each  z  €  0,  the  mappings  4  >-»  3{4i  z)  and  4  *-*•  Af(^,  z)  are  both  Frechet 
differentiable,  i.e.,  (A)  holds. 

Since  0  is  convex  and  the  mapping  K  is  linear,  we  have  that  for  any  4  €  z\  €  0, 

Z2  €  0  and  7  €  (0, 1),  the  element  z^  =  jzt  +  (1  -  7)22  belongs  to  0  and 

M(4,  z^)  =  A4  +  A  N(4)  +  A  Kz^ 

*  y(M  +  A  N(+))  +  (1  -  7)(*l#  +  A  *(#)) 

+  A(7ffzi  +  (1  -  i)Kzi)) 

=  7Af(^,xi)  +  (1  -  7)Af(^,z2)  • 

Moreover,  the  mapping  g  *-*  ||p||J  is  obviously  convex  from  0  into  R  so  that 

J(*,».)  =  |  #  -  *>118  +  |  111*  +  (1  ~  l)*48 

S  5  (7 II*  -  *48  +  (1  -  7)  II*-  *48)  +  ^(1  Ik  118  +  (1  -  7)  IMIS) 

Thus,  (B)  holds. 

The  operator  M^(^,g)  from  X  to  Y  is  defined  by 

M+{+,g)  +=A++  A  N'(4)  4  H§(U) , 

or  simply, 

Mi)(4,g)=A  +  \N'(4), 
where  for  any  4  €  X ,  the  operator  N'(4)  •  X  -*  Y  is  given  by 
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Thus,  using  the  definition  of  [•,•],  we  obtain  that  N'(if))  •  4  €  Ll(Cl)  =  Z.  Using  the  compact 
imbedding  Z  «-►«-»  Y,  we  see  that  JV'(^)  is  a  compact  operator  from  X  to  X.  As  a  result, 
=  A  +  A  =  A(I  +  A  A-1^'^))  is  a  Fredholm  operator  so  that  it  has  a  closed 

range  with  a  finite  codimension,  i.eM  (C)  holds. 

We  now  verify  that  the  hypothesis  (D)  holds  for  almost  all  A.  In  fact,  if  A  /  0  and  (1/A)  £ 
tr(-A-1  where  tr(-A-1  jV'(^))  denotes  the  spectrum  of  the  operator  (-A^N'W)),  then 

it  follows  that 

Y  =  Range{A(I  +  A  A~'N'(i>)))  =  Range{M^g)) 

so  that  Range(M^(i>, g))  contains  0  €  Y  as  an  interior  point,  i.e.,  (D)  holds.  Since  the  spectrum 
of  a  compact  operator  is  at  most  countable,  we  conclude  that  for  almost  all  A,  (D)  holds. 

Hence,  by  Theorem  3.1,  we  obtain  that  for  almost  all  A,  there  exists  a  if  €  Y*  =  Hq(SI)  such 

that 

A>i  +  A  [*'(*)]•  ■  t,  -  J^g)  =  0  (3.6) 

and 

A*,0,n)<  £(*.*.1)  Vr€0.  (3.7) 

Recalling  the  definition  of  the  operators  A  and  JV'(^)  as  well  as  the  fact  that  A  is  self-adjoint,  i.e., 
A*  =  A ,  we  see  that  (3.6)  is  equivalent  to  (3.1)-(3.2).  Using  the  definition  of  the  Lagrangian,  we 
easily  see  that  (3.7)  is  simply  a  restatement  of  (3.3),  the  latter  being  equivalent  to  (3.4).  Thus  we 
have  proved  that  a  Lagrange  multiplier  if  £  jfiT^(ft)  exists  such  that  (3.1)-(3.2)  and  (3.4)  hold.  □ 

Remark.  Now  we  can  explain  the  rescalings  that  were  introduced  in  §1.  With  these  rescalings, 
the  parameter  A  appears  in  the  form  of  an  eigen  value  in  the  operator  M^{ij},g)  =  A  +  A  N'(ip) 
so  that  spectrum  theories  for  compact  operators  can  be  readily  employed  in  the  proof  of  Theorem 
3.2.  0 

So  far  0  has  only  been  assumed  to  be  a  closed,  convex  subset  of  L2( fl).  In  the  sequel  we 
specialize  to  the  case  0  =  L2(tl).  Then  we  can  obtain  two  improvements  in  Theorem  3.2,  the  first 
being  that  the  results  of  Theorem  3.2  now  hold  for  all  A,  the  second  being  that  (3.4)  becomes  an 
equality. 

Theorem  3.3.  Assume  0  =  L2(fl)  and  £  Hq(H)  x  L2(S1)  is  an  optimal  solution  that 
minimizes  (1.8)  subject  to  (1.4)-(1.5).  Then,  for  all  A  €  A,  there  exists  a  Lagrange  multiplier 
i)  €  #o(Il)  satisfying  the  Euler  equations  (3.1)-(3.2)  and 

(t,g  +  rf j)  =  0  Vt€L2(fl).  (3.8) 

Proof.  In  view  of  Theorem  3.2  and  its  proof,  we  only  need  to  verify  (D)  when  (1/A)  6 
<r(-A-1JV'(u)),  i.e.,  we  need  to  show  that  for  each  /  £  H~2(ft),  there  exists  a  g  £  L2(il)  and  a 
i>  £  Bq(SI)  such  that 

A*  +  AJV'(*)  i>  +  \Kg  =  f, 

or  equivalently, 

o(&,*i)  +  \{[Th,H<h)  =  V<h€  HUH)  (3.9) 
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and 

a(i>%,<h)  -  -  A([0i,0a],02)  -  A (5,02) 

=  (/a,fc>  V<heH*(Sl). 


To  show  this,  we  first  let  0  €  H$(Sl)  be  a  solution  of 


(3.10) 


«(0i,0i)  +  A([02,02],0,)  =  (h,h)  V<fne  H*(Sl) 

and 

<*(02,02)-  A([0i,02],02)  =  (/a,  ^a)  V  0a  € 

The  existence  of  such  a  0  can  be  shown  in  a  manner  similar  to  that  for  showing  the  existence  of 
a  solution  to  the  von  Karman  equation;  the  key  step  is  that  by  adding  the  two  equations  with  the 
test  function  0  replaced  by  0,  we  have  the  a  priori  estimate 


<*(0i,0i)  +  <*(02,02)  =  </i,0i)  +  (/2,02>. 

Now,  having  chosen  such  a  0,  we  simply  set  g  =  — [0i,0a].  Note  that  regularity  results  for  the 
biharmonic  equation  applied  to  (1.4)-(1.5)  yield  0  €  H4((l)  (see  [3]).  Hence,  using  imbedding 
theorems  we  deduce  that  g  €  L2(Sl).  It  is  obvious  that  g  and  0  satisfy  (3.9)-(3.10).  Thus,  we  have 
verified  (D)  so  that  by  Theorem  3.1,  there  exists  a  i)  6  X  such  that  (3.1)-(3.2)  and  (3.4)  hold. 

Finally,  (3.4)  trivially  reduces  to  (3.8)  in  the  present  case  of  8  =  Z2(fl).  0 

Combining  the  results  of  §2  and  this  section  we  see  that  we  have  proved  the  existence  of  a 
triplet  (0,5,17)  €  J?o(^)  x  L2($l)  x  #oW  s**ch  that  the  system  (1.4)-(1.5),  (3.1)-(3.2)  and  (3.8) 
is  satisfied.  For  convenience  we  collect  these  equations  here  to  form  an  Optimality  System  of 
Equations : 

<*(0i,0i)  +  |((0J,02],0|)  =  0  V0i€#o(fi)> 

<*(02,02)  -  A([0i,02],0a)  =  X{g,<h)  V  02  €  Hq(SI), 

-  A ([02, 17a], Cl)  =  A(01  -  010, Cl)  VCi  €  J702(ft), 

a(C2,»?2)  +  A ([02, 171], C2)  -  A([0i,ij*],Ca) 

=  A(02  -  020, C2)  v  C2  e  f?o(ft) 
and 

(t,<7  +  i 72)  =  0  VteL*({l). 


(3.11) 

f  3.12) 

(3.13) 

(3.14) 

(3.15) 


4.  Finite  element  approximations  and  error  estimates 

4.1.  Definition  of  finite  element  approximations. 

A  finite  element  discretization  of  the  optimality  system  (3.11)-(3.15)  is  defined  in  the  usual 
manner.  For  simplicity,  we  will  only  study  conforming  finite  element  approximations  in  this  paper. 
However,  the  error  estimation  techniques  used  in  this  paper  are  equally  applicable  to  mixed  finite 
element  approximations  based  on  the  Hellan-Hermann-Johnson  scheme  for  biharmonic  equations. 
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See  [4]  (also  [5])  for  the  definition  and  discussions  of  the  mixed  Hellan-Hermann-Johnson  scheme 
for  biharmonic  equations. 

We  first  choose  families  of  finite  dimensional  subspaces  Xh  C  H q(Q)  and  Gh  C  Z»J(fl)  param¬ 
eterized  by  a  parameter  h  that  tends  to  zero  and  satisfying  the  following  approximation  properties: 
there  exists  a  constant  C  and  an  integer  r  such  that 

inf  II*  -  /||2  <  CAm||4||m+2  ,  V  4>  <E  ,  1  <  m  <  r  (4.1) 

4>k  exk 

and 

inf  \\z- zK\\o<Chm\\z\\m,  V  z  £  1  <  m  <  r.  (4.2) 

zh£Gk 

One  may  consult,  e.g.,  [2]  and  [6]  for  some  finite  element  spaces  satisfying  (4.1)  and  (4.2).  For 
example,  one  may  choose  Xh  =  Vh  x  Vh  where  Vh  is  the  piecewise  quintic-C1  (ft)  finite  element 
space  constrained  to  satisfy  the  given  boundary  conditions  and  defined  with  respect  to  a  family  of 
triangulations  of  ft.  In  this  case,  h  is  a  measure  of  the  grid  sue. 

Once  the  approximating  spaces  have  been  chosen,  we  may  formulate  the  approximate  problem 
for  the  optimality  system  (3.11)-(3.15):  seek  €  Xh,  gh  £  Gk,  and  rf1  £  Xh  such  that 


|  (« rtM)  =0  V  <t>x  £VK  ,  (4.3) 

~  A(K,^],^)  =  A  .(**,#)  V  $  6  V* ,  (4.4) 

<*«i\tf)  -  A([^,»?jfc],Cifc)  =  Ktf  -  V>i0,Cifc)  v  €  V* ,  (4.5) 

«(C2\rtf)  +  XWttflQ)  -  A([*\itf],tf)  =  V  €  Vh  (4.6) 

and 

(z\sA  +  J? 2)  =  0  Vz*€C?A.  (4.7) 


Note  that  in  the  last  equation,  if  Gh  =  Vh,  then  we  have  gh  —  -rjj  so  that  the  variable  gh  can  be 
eliminated  to  simplify  the  approximate  problem.  But  in  general  we  have  to  deal  with  the  entire 
system  (4.3)-(4.7). 

4.2.  Quotation  of  results  concerning  the  approximation  of  a  class  of  nonlinear  prob¬ 
lems. 

The  error  estimate  to  be  derived  in  Section  4.3  makes  use  of  results  of  [4]  and  [8]  concerning  the 
approximation  of  a  class  of  nonlinear  problems.  These  results  imply  that,  under  certain  hypotheses, 
the  error  of  approximation  of  solutions  of  certain  nonlinear  problems  is  basically  the  same  as  the 
error  of  approximation  of  solutions  of  related  linear  problems.  Here,  for  the  sake  of  completeness, 
we  will  state  the  relevant  results,  specialized  to  our  needs. 

The  nonlinear  problems  considered  in  [4]  and  [8]  are  of  the  following  type.  For  given  A  €  A, 
we  seek  ifr  £  X  such  that 

W(A,VO  =  V>  +  T0(A,^)  =  O,  (4.8) 

where  T  £  £(}>; X),  Q  is  a  C2  mapping  from  Ax  X  into  y ,  X  and  y  are  Banach  spaces,  and  A 
is  a  compact  interval  of  It.  We  say  that  {(A, V’(A))  :  A  £  A}  is  a  branch  of  solutions  of  (4.8)  if 
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A  — ►  0(A)  is  a  continuous  function  from  A  into  X  such  that  H(A,  0(A))  =  0.  The  branch  is  called 
a  regular  branch  if  we  also  have  that  A,  0(A))  is  an  isomorphism  from  X  into  X  for  all  A  G  A. 
Here,  •)  denotes  the  Frechet  derivative  of  H{-,  •)  with  respect  to  the  second  argument.  We 
assume  that  there  exists  another  Banach  space  Z ,  contained  in  y,  with  continuous  imbedding, 
such  that 

(?ifr(A,0)  €  C{X-,Z)  V  A  6  A  and  0  €  A* ,  (4.9) 

where  G^(',  •)  denotes  the  Frechet  derivative  of  £/(*,  •)  with  respect  to  the  second  argument. 

Approximations  are  defined  by  introducing  a  subspace  Xh  C  X  and  an  approximating  operator 
Th  G  £(y;«¥h).  Then,  given  A  €  A,  we  seek  0h  G  Xh  such  that 

=  +  =  (4.10) 

Concerning  the  operator  Tfc,  we  assume  the  approximation  properties 

lim  \\(Th  -  T)u\\x  =  0  wu>G  y  (4.11) 

/i— »o 

and 

lim\\(Tk-T)\\c(z.X)  =  Q.  (4.12) 

Note  that  whenever  the  imbedding  Z  C  ^  is  compact,  (4.12)  follows  from  (4.11)  and,  moreover, 

(4.9)  implies  that  the  operator  T<7^(A,  0)  €  £(/f;  A*)  is  compact. 

We  cam  now  state  the  result  of  [4]  or  [8]  that  will  be  used  in  the  sequel.  In  the  statement  of 
the  theorem,  D2G  represents  any  and  all  second  frechet  derivatives  of  Q. 

Theorem  4.1.  Let  X  and  y  be  Banach  spaces  and  A  a  compact  subset  of  K, .  Assume  that 
G  is  a  C2  mapping  from  A  x  X  into  y  and  that  D2G  is  bounded  on  all  bounded  sets  of  A  x  X. 
Assume  that  (4.9),  (4.11),  and  (4.12)  hold  and  that  {0(A);  A  G  A)  is  a  branch  of  regular  solutions 
of  { 4.8).  Then,  there  exists  a  neighborhood  O  of  the  origin  in  X  and,  for  h  <  ho  small  enough,  a 
unique  C2  function  A  — >  0h(A)  G  Xh  such  that  (0h(A);  A  G  A}  is  a  branch  of  regular  solutions  of 

(4.10)  and  0fc( A)  —  0(A)  G  O  for  all  A  G  A.  Moreover,  there  exists  a  constant  C  >  0,  independent 
of  h  and  A,  such  that 

!I0*(A)  —  0(A)||*  <  C^{Th  —  T)5(A,0(A))||*  VAGA.  □  (4.13) 

4.3.  Error  estimates. 

In  this  section  we  will  derive  error  estimates  for  the  finite  element  approximations  of  solutions 
of  the  optimality  system. 

We  begin  by  fitting  the  optimality  system  and  its  finite  element  approximations  into  the 
abstract  framework  in  §4.2. 

We  set  X  =  *  x  G  x  Y\  y  =  Y  x  X\  Z  =  Z  x  Z,  and  Xh  =  Xh  x  Gh  x  (Y*)fc.  Then 
we  have  that  Z  is  compactly  imbedded  into  y.  (We  recall  that  the  spaces  X  =  /fg(fl)  =  Y*, 
Y  =  H~3(Sl)  =  X *,  Z  =  Ll(l 2)  and  G  =  Z/2(fi)  were  all  introduced  in  §3.)  We  define  T  G  C(y;X) 
a 8  follows:  T(f,f)  =  (0,  g,  7j)  for  (f,f)  G  and  (0,j, »))  G  A”  if  and  only  if 

a(0,0)  =  (f,0)  V  0  G  ^03(n) ,  (4.14) 
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o(^,w)  =  (f,u)  Vue  H$(Sl) 


(4.15) 


and 

(*,g  +  m)  =  o  vteL*(ii).  (4.16) 

Similarly,  the  operator  Th  e  C(y;Xh)  is  defined  as  follows:  Th(f,f)  =  (^^ff4,*)4)  for  (r,r)  G  y 
and  (0A,fffc,ijk)  6  «¥*  if  and  only  if 


0(0  V)  =  (f,0*) 

V 0A  €  Af *, 

(4.17) 

o(ihM)  =  (r,vk) 

v^ex* 

(4.18) 

(th,gh  +  Tft)  =  0 

vtkeG*. 

(4.19) 

The  nonlinear  mapping  Q  :  A  x  X  — ►  y  is  defined  as  follows:  Q (A,  (0,0,i|))  =  (r,r)  for  A  €  A, 
(0>9>il)  €  A’,  and  (f,f)  €  y  if  and  only  if 


(  -*(02,  &]  \ 

>&)  +  *#/ 


and 

z.  _  (  _  +  A(0i  -  0m)  ^ 

\-A(ife,^]  +  A[0i,tfi]  + A(02 -020)/  ’ 

or  equivalently, 


{c(K(+,s,i)),<+,<)) 

=  -^([Vo,Vol,4»i)  +  +  HSi'h) 

+  A([02,ija],Ci)  +  A(0i  -  0io,<i)  -  A($3,i7i],Ca) 

+  A([0i,i»],<2)  +  A(0,  -  0ao, Ca)  V  (0,0  e  BZ(Sl)  x  F0a(fl). 


Using  the  operators  N  and  K  introduced  in  §3,  we  can  also  write  the  defining  equations  for 
Q  (A,  (0,  g,  if))  as  follows: 


f  =  \N(1>)  +  \Kg 
and 


(4.20) 


f  =  A[AT'(0)]*^-A(0-0o). 


(4.21) 


It  is  easy  to  see  that  the  optimality  system  (3.11)-(3.15)  and  its  finite  dimensional  counterpart 
(4.3)-(4.7)  can  be  written  as 

(1>,g,vi)  +  TS{\,{+>g>v))  =0 


and 

(*\  »*,!)*)  +  r^A^V,^))  =  o, 

respectively,  i.e.,  with  ip  —  (th,p,t l)  and  ipk  =  (iph,gh,Vh)>  in  the  form  of  (4.8)  and  (4.10),  respec- 
tivdy. 
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We  now  turn  to  the  verification  of  the  assumptions  in  Theorem  4.1.  We  first  examine  the 
approximation  results  for  the  linear  operators  T  and  Tk. 

Lemma  4.2.  For  every  (f,f)  €  y, 

/m||(T- T*)(f, r)||*  =  0. 

h- *0 

Furthermore,  ifT(f,f)  =  (^,#,9)  €  Pm+2(H)  x  Hm{Sl)  x  Hm+2(Sl),  then 

II (T  "  Tk)(f, f)|U  <  Ch*n(||^||m+2  +  IlffHm  +  MU+S)  • 

Proo/.  For  each  given  (f,f)  g  let  d>,g,n)  =  T(f,f)  and  ,gh,nh)  =  Th{f,f),  i.e., 
€  X  and  ,gk,^h)  €  satisfy  (4.14)-(4.16)  and  (4.17)-(4.19),  respectively.  From  the 
well-known  results  concerning  the  approximation  of  biharmonic  equations  (see,  e.g.,  [2]  and  [6]), 
we  obtain  that 

lim  ||V>  -  ^||a  =  0 

/i — e-0 

and 

(tmP-rt  =  0.  (4.22) 

Furthermore,  if  ^  g  Hm+3(tl)  and  ^  g  Hm+2(Sl),  then 

%*-+kb<ch»mm+i 


and 

II#  -  V%  <  Cfcm||#||m+2  • 

(4.16)  and  (4.19)  yields  ( g  ~  gh,zh)  =  -(%  —  rfi,  zh)  for  all  zh  g  Gh.  We  then  have  that 

ll$-$hllo  =  (9-  gh,g-gh) 

=  (g  -  gh,g  -  zh)  +  (g  -  gk,zh  -  gh) 

=  (g  -  gh,g  -  *h)  -  (%  -  v* .  z1*  ~  9h ) 

£  II#  -  **118  +  jllj  -  S*llo  +  jllft  -  %‘llo 

+  |ll*‘  -  Silo  +  llfr  -  iillo  +  jIIj  -  S*llo 

so  that 

11$  -  flfc||o  <  c  (||S  -  ar^llo  +  11%  -  V2  ||o)  • 

It  follows  from  (4.2)  and  (4.22)  that  Iitoa-mj  ||jr  -  =  0.  Hence, 

urn  hc r  -  r»)(  f,f)\\x  =  u m  (y*  -  /y,  +  II  j  -  s‘||o  + ||«  -  «*|,)  =  o . 

Furthermore,  if  (f,g,ti)  g  27”*+2(Q)  x  Hm(Sl)  x  J5Tm+2(fi),  then 

IIS -  SNIo  <  Chr"(\\g\\m  +  ||%|w+2) 
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so  that 


-A. 


and 


II (T  -  T‘)(f,r)|U  =  ||*  -  *-||,  +  ||j  -  5*110  +  ||i»  -  $#||i 

<  Cfc"-(|W|„M  +  ||j|U  +  NU+j)-  0 

Next,  we  examine  the  derivatives  of  the  mapping  Q. 

Lemma  4.3.  The  mapping  G  is  twice  continuously  Prechet  differentiable;  the  second  order 
derivative  is  bounded  on  all  bounded  subsets  of  Ax  X;  and  the  first  order  derivative  maps  X  into 
2  continuously,  i.e.,  €  C(X\Z). 

Proof.  Using  imbedding  theorems,  Holder  inequalities  and  the  fact  that  G  consists  of  poly¬ 
nomial  maps,  we  deduce  that  the  mapping  G  is  twice  continuously  Frechet  differentiable  and  the 
second  order  derivative  is  bounded  on  all  bounded  subsets  of  A  x  X. 

A  simple  calculation  shows  that  g  ^  €  £(X;y)  is  given  by 

Grf,g,il)  (A’  9 1  9'®  =  A  ( i  -  i> ) 

where  the  operators  N  and  K  along  with  their  derivatives  were  introduced  in  §3.  Sobolev  imbedding 
theorems  imply  that 

€!•(«)■ 

Of  course,  1>  6  L1(ft).  From  the  definition  of  the  operator  A  we  see  that  K  maps  L2( ft)  into 
{0}  x  L2{Sl)  C  Ll(ft).  Hence  the  operator  ^  (A,  (^,  </,»?))  maps  X  into  2.  Moreover,  using 
Holder  inequalities  and  Sobolev  imbedding  theorems  it  is  easy  to  see  that 

l|G(5[#.,p(M^>J»^)) '  (^>S.$)IU><n> 

<  C(|M>  +  Wli  + 1)  (M>  +  ||5l|o  +  ») 

so  that 

Next,  we  recall  the  notion  of  regular  solutions.  A  solution  of  the  optimality 

system  (3.11)-(3.15)  is  called  regular  if 

(/  +  T G(4>,  9,*l)(*A4>,9,  n)))  •  (&  9,1)  =  (l  k  0 

has  a  unique  solution  6  X  for  any  ((,  h,Q  €  X.  As  T  is  invertible,  ((,  h,Q  =  T(f,f)  with 

(f, f)  €  y  =  x  Thus,  (^(A),jf(A),iy(A))  is  a  regular  solution  of  the  optimality 

system  (3.1 1)-(3.15)  if 

(/  +  =  T(f,f) 
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has  a  unique  solution  (^,£,9)  €  X  for  any  (f,f)  £  y,  or  equivalently,  the  following  system  (for 
the  unknowns  (^,  <j,ty))  is  uniquely  solvable  for  any  (f,f)  £  y-. 


a(V>i,<fo)  +  1)  =  (n,<fa)  V <fn  G 

(4.23) 

a{rfa,fa)-  A([i(n,rh\,<h)  -  A(tyi,t&2},<fc)  -  A(g,<h) 

=  (ri,<M  V<fcef?02(fi), 

(4.24) 

°((j»»h)  -  M(V>3>$»],Ci)  ~  -  MV>i*Ci) 

=  (n,Ci>  VCi€fl0J(ft), 

(4.25) 

a(Ca,th)  +  A([V>2, ^1], Ca)  +  ^(['02, »7i], C2)  -  ^ ([V'l , ^2], C2) 

-  A([t^i,i72],Ca)  -  KM)  =  (fa,  (2)  v  C2  €  /f02(«) 

(4.26) 

and 

M  +  »h)  =  0  Vt€L2(n). 

(4.27) 

Note  that  the  linear  operator  appearing  on  the  left  hand  side  of  (4.23)-(4.27)  is  obtained  by 
linearizing  the  optimality  system  (3.11)-(3.15)  about  (ip,g,ri). 

Lemma  4.4.  For  almost  all  A,  solutions  (^(A), 5(A),  17(A))  of  the  optimality  system  (3.11)- 
(3.15)  are  regular. 

Proof.  The  system  (4.23)-(4.27)  can  be  rewritten  as 

(/  +  ATS(1>,g,rlj)(i>,g,rl)  =  T(f,f) , 

(4.28) 

where  the  linear  operator  S(1>,g,rj)  :  X  -*  y  is  defined  by 

?  N'W-i+Kg  \ 

It  was  established  in  Lemma  4.3  that  the  mapping  Q ^  g  ^  (A,  (^,47,1}))  is  compact  from  X  into  y. 
Now,  T  €  C(y\X),  so  that  (/  +  A  TS(rf>,g,ri))  is  a  compact  perturbation  of  the  identity  operator 
from  X  to  X.  Thus,  for  almost  all  A,  (4.28),  or  equivalently  (4.23)-(4.27),  is  uniquely  solvable,  i.e., 
for  almost  all  A,  the  solution  (tf(A),jr(A),  17(A))  of  the  optimality  system  (3.1 1)-(3.15)  is  regular. 
0 

We  are  now  in  a  position  to  derive  error  estimates.  In  the  following  theorem,  we  will  assume 
that  the  solution  (^(A),j(A),  17(A))  of  the  optimality  system  (3.11)-(3.15)  is  regular.  Lemma  4.4 
guarantees  that  this  is  indeed  the  case  for  almost  all  A.  Lemmas  4.2  and  4.3  verified  all  the 
assumptions  in  Theorem  4.1.  Thus  we  are  led  to  the  following  error  estimates. 

Theorem  4.5.  Assume  that  A  is  a  compact  interval  0/  R  +  and  that  there  exists  a  branch 
{(^( A), 17(A),  17(A))  :  A  €  A}  of  regular  solutions  of  the  optimality  system  (3.11)-(3.15).  Assume 
that  the  finite  element  spaces  Xh  and  Gh  satisfy  the  hypotheses  (4.1)-(4.2).  Then,  there  exists  a 
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6  >  0  and  an  h o  >  0  such  that  for  k  <  ho,  the  discrete  optimality  system  (4.3)-(4.7)  has  a  unique 
branch  of  regular  solutions  {{^h(X),gk(X),tfh(X))  :  A  e  A}  satisfying 

{||**(A)  -  *A)||a  +  ||sffc(A)  -  ?(A)||o  +  ||^(A)  -  *KA)||2}  <6  for  all  A  €  A. 

Moreover, 

lim{ ||*fc(A)  -  «A)||,  +  ||^(A)  -  ff(A)||o  +  ||^(A)  -  r^X)\\2}  =  0 , 
uniformly  in  A  €  A. 

If,  in  addition,  the  solution  of  the  optimality  system  satisfies  (i(t{X),g(X),r){X))  6  Hm+2(Sl)  x 
x  Hm+2(fl)  for  A  €  A,  then  there  exists  a  constant  C,  independent  of  h,  such  that 

INK  A)  -  *‘(A)||,  +  ||,(A)  -  j»(A)||o  +  WA)  -  t|A(A)||j 
<  C*»(WA)S«+,  +  ||s(A)||m  +  |fo(A)||m+,) , 

uniformly  in  A  €  A.  0 


5.  Conclusions 

We  studied  an  optimal  control  problem  for  the  von  Karman  equations  in  this  paper.  We  first 
gave  the  mathematical  statement  of  the  problem  and  proved  the  existence  of  an  optimal  solution. 
We  then  applied  the  Lagrange  multiplier  rules  to  derive  an  optimality  system  of  equations  that  the 
optimal  state  must  satisfy  (the  use  of  Lagrange  multiplier  rules  was  justifed).  We  finally  defined 
finite  element  approximations  for  the  optimality  system  and  derived  optimal  error  estimates.  The 
functional  we  minimized  was  a  tracking  functional  (i.e.,  the  tracking  of  the  Airy  stress  function 
and  the  normal  deflection).  The  control  we  used  was  a  distributed  control  (i.e.,  the  external  load). 
However,  the  methods  used  in  this  paper  apply  equally  well  to  optimal  control  problems  with  other 
objectives  (e.g.,  minimizing  the  stress  functional  in  some  areas)  and/or  other  types  of  controls  (e.g, 
boundary  controls). 
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